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Infrared  Chemiluminescence  from  Oxygen 
Atom  Attack  on  Fuel  Species. 

R  H.Kummler  and  B  R  Fisher 


INTRODUCTION 


The  energy  released  in  an  exothermic  chemical  reaction  can 
take  several  forms  depending  on  the  reaction  path.  The  chemical 
energy  of  the  reactants  can  be  converted  into  translational 
energy  of  the  products  or  it  can  be  converted  into  internal 
energy  of  the  products  or  some  combination  of  both.  The 
internal  energy  can  be  distributed  among  the  electronic, 
vibrational,  or  rotational  modes  of  the  product  molecule  or 
molecules.  If  the  product  is  elevated  to  a  quantum  state  which 
has  an  optically  allowed  transition  to  a  lower  state,  then  the 
resulting  light  emitted  permits  the  study  of  the  reaction  path 
leading  to  the  internal  energy.  The  observation  of  radiation 
associated  with  this  non-equilibrium  process  thus  provides 
information  regarding  the  fundamental  chemical  kinetics  of 
the  reaction  (e.g.,  Charters  and  Polanyi  (1960);  Polanyi,  1966: 

Pacey  arid  Polanyi,  1971;  Kushafar,  et  al,  1971).  It  also 
provides  the  opportunity  to  assess  the  presence  of  one  or  both 
of  the  reactant  partners  in  an  unknown  mixture  of  reactant  gases. 

This  capability  is  of  considerable  interest  in  the  air  pollution 
monitoring  field  and  is  applicable  to  both  local  and  remote  monitoring 
of  exhaust  gases  (of  rockets  as  well  as  automobiles)  in  the  presence 
of  an  oxygen  atom  enviroment.  In  the  spectral  region  of  wave¬ 
lengths  from  ljji  to  roughly  20  or  30 jx,  the  emitted  light  generally 
characterizes  vibrational  transitions  which  can  result  when  the 
exothermic  reaction  produces  a  new  bond  in  a  heteronuclear  molecule. 


In  previous  work  (Krieger,  Malki,  and  Kummler,  1972;  Malki,  1972; 
Krieger  and  Kummler,  1973)  we  have  investigated  the  light  emitted  by 
reactions  which  might  be  expected  to  occur  under  conditions  appropriate 
to  hydrocarbon  combustion  in  the  upper  atmosphere.  First,  we 
identified  the  potential  radiatiors  on  the  basis  of  the  known  ground 
state  chemistry  and  energetics.  Second,  we  identified  the  visible 
and  near  infrared  (7000-9000A)  emission  which  illustrated  the 
presence  of  both  electronically  and  highly  vibrationally  excited  OH 
early  in  the  reactor  time  scale.  Third,  we  estimated  the  quantum 
yield  for  the  emission  of  photons  at  2.7y  (Kummler,  Fisher,  and 
Malki,  1973)  using  non-dispersive  techniques  for  0+C2H4.  While 
the  indirect  evidence  appeared  to  strongly  support  OH  as  the  radiating 
species  in  all  of  the  previous  work,  we  had  not  been  able  to  obtain  the 
direct  spectral  evidence  needed  to  positively  identify  the  radiator. 

In  this  work,  we  have  concentrated  on  providing  the  spectroscopic 
measurements  necessary  to  achieve  the  identification  for  a  variety 
of  potentially  important  systems: 


(1) 

o+c2h4 

(6) 

0+C2H6 

(2) 

0+c2h2 

(7) 

0+C4H10 

(3) 

h+no2 

(8) 

0+C3H6 

(4) 

O+HCHO 

(9) 

OH+CO 

(5) 

h+ho2 

(10) 

oh+c2h4 
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EXPERIMENTAL  PROCEDURE 


All  experiments  were  conducted  at  roughly  1  Torr  total 
pressure  using  a  discharge  flow  system  shown  schematically  in 
Figure  (1).  This  system  is  very  similar  to  that  previously 
decribed  by  Krieger,  Malki  and  Kummler  (1972)  and  by  Malki  (1972). 

For  convenience  we  will  repeat  that  description  as  follows. 

The  flow  tube  reactor  is  a  4  ft.  long  stainless -steel  tube 
with  a  1  inch  I.D.  teflon  tube  liner  which  aids  in  preventing 
wall  recombination  of  some  atomic  species.  It  is  provided  with 
six  ports  for  radial  viewing  through  quartz  windows  (Muffaletto 
Optical  Company,  Baltimore,  Md.)  Two  additional  ports  were  used 
for  measuring  the  absolute  pressure  and  the  drop  in  pressure  across 
the  reactor.  Linear  average  velocities  up  to  80m/sec(with  70  m/sec 
being  a  typical  velocity)  were  achieved  when  the  pressure  in  the 
reactor  was  about  1  torr  using  a  Heraeus  Englehard  (E - 2 2 5  ,  147  cfm) 
mechanical  pump  and  Roots  Blower. 

A  Baratron  (MKS  type  77)  capacitance  manometer  with  a  three  Torr 
head  was  employed  to  monitor  the  reactor  pressure  to  better  than  O.ly 
accuracy.  An  alphatron  (type  520)  and  a  transducer  (Viatran  model 
304,  0-15  psia)  were  the  pressure  sensing  devices  used  to  monitor 
ballast  tank  pressures.  Their  respective  locations  and  functions  in 
the  system  are  illustrated  (as  well  as  other  equipment)  in  Table  (1). 

The  detectors  used  for  this  work  are  the  PbS  detector  (Santa 
Barbara  Research  ITO  type)  in  conjunction  with  a  phase  sensitive  lock-in 
amplifier  (PAR  HR-8)  for  IR  detection,  and  the  photomultiplier  detector 
(RCA  C31025C)  in  conjunction  with  an  electrometer  (Keithley 602)  for 
detection  of  visible  emission. 
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In  order  to  get  absolute  emission  values,  the  system  had  to 
he  calibrated  using  the  reaction  (Fontijn  et  al,  1964) 

NO  +  0  N02  +  hv  (2.7 y) 

The  actinometric  standard  of  Fontijn,  et  al  was  extended  into  the 
2.7y  region  using  the  data  of  Kaufman  (1973). 

The  oxygen  atoms  were  generated  by  passing  oxygen  (Cryogenic 
sales  99.6%)  pure  or  diluted  in  He  (Cryogenic  sales  99%)  through 

an  Evansen  type  cavity  discharge  powered  by  a  2450  MHz  Raytheon 
PGM  10X2  diathermy  unit. 

The  oxygen  atoms  were  passed  into  the  flow  tube  reactor  through 
a  closed  end  teflon  tube  with  a  radially  multi-perforated  exit 
stream.  This  permitted  thorough  mixing  and  a  plug  flow  pattern  at 
the  point  of  nixing  due  to  the  turbulence  caused  by  introducing  the 
gas  through  a  "showerhead"  gas  inlet.  The  plug  on  the  end  of  the 
teflon  piece  served  also  to  prevent  the  discharge  region  radiation 
from  reaching  the  observation  ports  through  scattering.  The  oxygen 
atom  concentrations  were  measured  by  using  the  NC>2  titration 
technique  (Kaufman  1964  ). 

A  steady  state  flow  of  the  gas  typical  of  the  experiment  (02  or 
02  *  He  or  H2  *  He)  was  obtained.  Then  the  microwave  generator  was 
turned  on.  Argon  or  Helium  is  usually  used  to  keep  the  free  atoms 
and  radicals  at  low  combinations  ,  thus  preventing  recombination. 

It  helps  stabilize  the  discharge  and  it  also  served  to  keep  the  pressure 
in  the  reactor  constant  allowing  a  reasonably  constant  linear  velocity 
upon  adding  the  ground  state  gas. 
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Helium  contributes  most  of  the  pressure  in  the  reactor 
Oi  torr)  while  the  ground  state  stream  (NO,  N02 ,  C2H4,  etc.)  and 
the  02  or  H2  are  relatively  small  fractions  of  the  total  pressure 
The  phase  angle  on  the  lock-in  amplifier  changed  as  one 
moved  the  PbS  detector  downstream  and  as  the  pressure  changed  in 
the  flow  tube.  This  had  to  be  checked  for  all  experimental  data 
points.  The  reference  to  the  PAR  was  taken  from  the  pulse  generator 
of  the  PGM  10X2.  The  pulse  rate  of  the  discharge  was  adjusted  so 
that  a  single  pulse  of  atoms  traveled  down  the  flow  tube.  This 
typically  required  repetition  rates  of  34  Hz  to  100  Hz.  The  duty 
cycle  for  most  measurements  was  501. 

A  pressure  base  line  and  a  zero  intensity  base  line  were  recorded 
on  a  two  channel  Honeywell  electronic  ? 94  recorder.  As  the  ground 
state  gas  was  added  to  the  reactor,  the  intensity  of  the  light  emitted 
from  the  excited  species  was  recorded  with  the  corresponding  preset 
scale  on  the  lock-in  amplifier.  Utilization  of  the  capillary  flow 
meter  permitted  direct  measurements  of  the  intensity  vs.  upstream 
pressure  at  the  capillary  using  the  PAR,  the  Viatran  pressure  trans¬ 
ducer  and  a  Honeywell  Model  (56))  XY  recorder. 

A  run  that  involved  the  reaction  NO  +  0  ■+■  02  +  hv  accompanied 
every  set  of  data  to  put  tL '  lock-in  amplifier  readings  on  an 
absolute  basis.  A  titration  for  oxygen  atoms  was  also  done  as  a 
part  of  the  system  calibration.  The  PbS  detector  with  the  PAR  HR8 
was  used  to  monitor  the  intensity  peak  of  the  reaction  NO,  *  0  ■+•  NO 

M 

+  02  as  a  function  of  N02  concentration. 
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THE  MECHANISMS 


A.  0  +  C2H4  and  0  «■  HCHO 

The  mechanism  for  oxygen  atom  attack  on  ethylene  is  well 
established  in  terms  of  ground  state  chemistry  and  Malki  (1972) 
has  investigated  the  emission  at  2 . 7 p  in  considerable  detail. 

The  reaction  of  oxygen  atoms  with  ethylene  has  been  phenomeno¬ 
logically  found  to  produce  infrared  radiation,  at  least  partly 
in  the  2 . 7  p  region  (Malki,  1972).  The  major  effort  of  this  work  is 
devoted  to  the  qualitative  explanation  of  that  radiation  in  terms 
of  detailed  spectra  on  which  to  base  molecular  mechanisms.  The 
general  ground  state  system  of  intermediates  resulting  from  the 
combination  of  0  and  C2H4  is  well  established  primarily  from  the 
discharge-flow-tube-time-of -flight  mass  spectrometric  study  of  Niki 
et  al  (1969)  and  the  direct  observation  of  the  products  CH3  and  CHO 
m  the  high  intensity  crossed  beam  experiment  of  Kanofsky  et  al  (1972) 
These  two  studies  along  with  the  similar  flow  tube  mass  spectrometric 
work  of  Herron  and  Penzhorn  (1969)  have  established  that  the  primary 
elementary  step  (901)  in  the  reaction  mechanism  is 

0  +  C2H4  -*■  CHj  ♦  CHO  (1) 

The  eventual  stable  products  formed  by  secondary  reactions  include 
H2C0,  CO,  H2,  02,  and  H.  Consideration  of  the  temporal  behavior 
of  these  products  as  well  as  isotopic  experiments  led  Niki,  et  al 
(1969)  to  postulate  that  the  following  reactions  occur: 

CHj  +  0  HCHO  +  H 
CHO  +  0  •♦•CO  +  OH 
HCHO  +  0  -*■  CHO  ♦  OH 
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(2) 

(4) 

(17) 


CHO  +  0>C02  +  H 

(5) 

OH  +  0  -*■  02  +  H 

(6) 

CHO  +  H-*-H2  +  CO 

(7) 

Water  was  not  observed  to  be  a  product  of  this  system  by  Niki 
but  his  system  employed  very  small  concentrations  of  C2H^ 

0  in  excess,  so  that  the  OH  concentration  and  hence  the  rate  of  the 
reaction: 

OH  +  oh-*-h2o  +  0 

were  both  small.  In  a  more  general  case,  H20  might  be  a  significant 
product,  and  it  could  be  a  significant  radiator  as  well. 

Both  OH  and  H20  are  produced  by  exothermic  reactions  in  the 
above  sequence  under  conditions  where  each  is  the  new  bond  formed 
by  a  neutral  rearrangement  reaction.  Under  conditions  such  as  this, 
it  is  well  established  that  a  significant  fraction  of  the  exothermicity 
of  the  reaction  can  appear  as  vibrationa.i  energy  in  the  new  bond 
(Smith  (1972);  Pacey  and  Polanyi  (1971);  McGrath  and  Norrich  (1957), 
Parker  and  Pimentel  (1969);  Jonathan,  et  al  (1971);  Stair  and 
Kennealy  (1967);  Hushfar,  et  al  (1971)  and  many  others). 

The  Basic  mechanism  is  illustrated  in  Figure  2  and  detailed 
in  Table  1 1 . 
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TABLE  II  (from  Malki,  1972) 


Reaction 

Rate 

Constant ; 
300°K 

Reference  AH 

1.  C2H4  ♦  O  +  CHj  ♦  CHO 

8 . 42 ' 12 

Davis,  et  al 

2.  CHj  +  0  -*■  HCHO  +  H 

3. O'11 

(1972) 

Niki,  et  al 

3.  HCHO  +  O+CHO  +  OHf 

(1969) 

Calculated 

4.  CHO  +  O  +  OH 

1.0‘12 

N’ki,  est. 

5.  CHO  +  0  C02  ♦  H 

3. O'13 

Niki,  est. 

6.  OH  ♦  0-02  ♦  H 

5 . 011 

Kaufman (1964) 

7.  CHO  ♦  H-  CO  ♦  H_ 

rsi 

H 

l 

O 

• 

Niki,  ct  al 

8.  OHf(v.1}-OH  *hv 

9.  CHO  ♦  02  CO  +  H02 

3.3 

(1969) 

Potter,  et  al 

2X10’13 

(1971) 

10.  OHf+  0  +02  +  H 

5. O'11 

Kaufman  (1964) 

11.  OH  +  OH  -  H20  ♦  0 

2 . 5  12 

Pixon-Lcwis  et 

12.  OH  ♦  OH  H2Of-  ♦  0 

2 . 5  13 

al  (1966) 

Estimate 

13.  CH^  ♦  02  ♦  M •♦  CHj02  ♦  M 

8. O'32 

Heicklcn 

14.  H  ♦  02  ♦  M-*  H02  ♦  M 

3. O'32 

Kaufman  (1964) 

15.  H2Of  ♦  M-H20  ♦  M 

1— • 

• 

o 

• 

1— • 
NJ 

Estimate 

16.  H2°  *H2°+hV 

30 

Pcnncr  (1959) 

17.  HCHO  ♦  0  •♦CHO  ♦  OH 

Niki,  ct  al 

18.  CHO  ♦  0-  CO  ♦  OHf* 

(1969) 

Calculated 

19.  OH  ♦  HCHO-H2Of+  CHO 

1.4'12 

Estimate 

20.  OH  ♦  H2  ♦  H .0  ♦  H 

6.6’1S 

Greiner (1969) 

k  cal 
gmole 

42.45 

66.54 

73.09 

98.05 

-16.79 

-75.62 

-18.88 

-16.95 

-47.4 

-24.77 


Estimate 

Greiner  (1969)  -119.234 


21. 

OH  +  H02  H20  +  02 

i.o-11 

Kaufman  (1964) 

-71.87 

22. 

H02  +  0-*-0H  +  02 

1.0'11 

Lloyd  (1970) 

-54.2 

23. 

H  +  H02  *►  OHf  +  OH 

I.O'12 

Crutzen  (1971) 

-38.4 

24. 

OH  +  H02  -*•  H20  f+  02 

l.O’12 

Estimate 

25. 

C2H4  +  0  -  CH3  +  CHOf 

5 . 2X1014 

Estimate 

26. 

CHOf  +  0  -►  CO  +  OHf 

l.O-11 

Estimate 

27. 

H  +  HO  -►  OH  +  OH 

l.O'11 

Estimate 

28. 

0  +  HOz-*-OHf  +  02 

l.O'12 

Estimates 

29. 

CO  +  0H-C02  +  H 

1.5-13 

Wilson  (1972) 

30. 

OHf  +  02  -02(b1Eg)  OH 

3X1014 

Potter,  et  al 

31. 

Ciff  +  H  -*•  H2  +  0 

1X1. O"12 

Estimate 

32. 

°2  (bl  E)  *02-02.02 

3X1016 

Estimate 

33. 

C2  (bl  +  °-*°2  *  0 

2X10"13 

Estimate 

34. 

02  (b1  E)  +  H-02  +  H 

2X1013 

Estimate 

35. 

H  +  HCHO  -►  H2  +  CHO 

1X10’13 

Westenburg  and 

36. 

CHOf  +  02  -  CO  +  H02 

1X10"11 

Dehaas  (1972) 

Estimate 

37. 

0-1/2  02 

137 

Estimate,  wall 

38. 

H-l/2  H2 

137 

deactivation 

Estimate 

39. 

OHf  -v  OH 

137 

Estimate 

40. 

H£0f-H20 

137 

Estimate 

41. 

CHOf  -  CHO 

137 

Estimate 

42. 

02  (b1  Eg)  -J2  (X3E) 

137 

Estimate 

43. 

44. 

Dummy 

OHf  +  M-OH  +  M 

1.5X10’14 

Estimate 

45. 

OH  f  +  H2  H20  +  H 

1X10‘14 

Estimate 

46. 

C2H4  +  OH-CHj  +  HCHO 

1 . 8X10' 12 

Morris  6  Niki  | 

Cl  971 ) 

In  the  general  case,  as  in  Figure  2,  both  the  OH  fundamental  and  the 
H2O  Vj  band  could  be  expected  to  contribute  to  emission  near  2.7p. 

In  order  to  simplify  the  initial  kinetic  system,  Malki  (1972) 
decided  to  work  at  small  C2H4  concentration  in  order  to  minimize 
the  H2O  production.  In  addition,  the  02  concentration  was 


minimized  in  order  to  avoid  the  reactions 

H  +  02  +  M-*-HC2  ♦  M  (14) 

0  +  H02  -*-OH+  +  02  (29) 

H  +  H02  -0H+  ♦  OH  (24) 

OH  +  H02  -  H2Of  ♦  02  (25) 

The  minimization  of  the  also  prevented  significant  contributions 

from 

OH  +  c2H4  products  (47) 

OH  +  H£  -  H20+  ♦  H, 

as  well  as  minimizimg  tnc  contribution  of 

0H+  +  C2H4  products  (48) 

0H+  ♦  HC  -  fragment  products 


to  the  OH  loss.  Under  such  circumstances,  the  major 0H  loss 
process  will  be  via 

0H+  +  0  -►  02  >  H,  (10) 

which  is  extremely  fast,  even  for  the  ground  state  OH  molecule. 

In  fact,  reaction  with  0  as  a  homogeneous  loss  process  is  so 
rapid  that  OH + deact ivation  on  the  wall,  which  would  be  diffusion 
controlled  even  if  every  wall  collision  were  effective  in  de¬ 
activation,  can  be  neglected.  In  the  present  work  the  entire 
range  of  stoichiometric  possibilities  were  covered  in  order  to 
test  all  of  the  above  hypothesis. 
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B.  H  ♦  N02 

The  reaction  products  for  this  system  have  been  established 
by  Kaufman  to  be  NO  ♦  OH.  Energetically,  only  two  quanta  of 
vibration  can  be  given  to  OH  in  the  creation  of  that  new  bond. 
Kaufnun  accounted  for  nearl>  all  of  the  product  in  the  ground 
vibrational  <=tate  through  absorption  spectroscopy.  We  had 
previously  investigated  this  reaction  via  emission  spectroscopy 
and  agreed  that  the  quantum  yield  for  vibrationally  excited 
OH  was  very  small.  A  complicating  factor  is  the  reaction 
OH  ♦  OH  -  H20  ♦  0, 

in  which  the  product  H20  could  also  radiate  in  the  2.7u  region. 
Hence,  radiometric  measurements  alone  at  2.7p  could  not  distinguish 
between  OH  and  H.,0  emission,  although  kinetic  considerations  would 
favor  OH.  The  rate  constant  for 
H  ♦  N02  -  NO  ♦  OH 

is  4.9  X  10*11  cc/sec  at  300°K,  whereas  the  rate  constant  for 
OH  ♦  OH  H20  ♦  0 
is  2.5  X  10  ^  cc/sec. 
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C.  0  +  C2H2 


The  product  of  oxygen  atom  attack  on  acetylene  has  been 
established  by  Kanofsky,  et  al  (1972)  and  by  Jones  and  Bayes  (1972) 
to  be 


o  +  c2h2  -*-ch2co 


ch2  +  CO 
HCCO  +  H 


c20  ♦  H2 


251 

15  -  4 5 %  l  751 
30  -  601 
.3% 

with  HCCO  being  rapidly  consumed  by  0  atoms  to  give 
HCCO  +  0-H+2C0 


Thus,  CO  is  a  major  species  formed  by  very  exothermic  reactions 
so  that  highly  vibrationally  excited  CO  is  expected  from  his 
system,  but  little  or  no  OH  or  H20  is  expected. 


D.  H  +  H02 

The  reaction 

H  ♦  H02  ■*>  20H+ 

has  been  studied  by  Charters  and  Pnlanyi,  (1960).  Those  authors 
concluded  that  this  very  rapid  reaction  produces  vibrationally 
excited  OH't 
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INFRARED  SPECTOSCOPY:  RESULTS  AND  DISCU S S I ON 

In  rder  to  establish  the  mechanism  for  0  atom  attack  on 
reactive  hydrocarbons  and  thereby  quantify  and  interpret  filter- 
kinetic  data,  in  both  the  visible  and  infrared,  we  have  taken 
extensive  spectral  measurements  in  the  1-3U  region,  using  a 
Jarrel  Ash  quarter  meter  monochrometer  and  an  ITO  Pbs  detector. 

The  total  pressure  for  these  measurements  was  about  1000-  1200y. 

The  wavelength  calibration  for  this  region  was  achieved 
using  the  CO  first  overtone  generated  by  the  reaction  of  0  +  C  H 
illustrated  in  Figure  3.  For  most  of  the  spectra,  the  resolution 
was  better  than  0.03^  so  that  the  location  of  the  maxima  between 
2.33  and  2.61  provides  an  excellent  calibration. 

For  the  oxygen  atom  attack  on  ethylene,  three  prominent  features 
are  observed.  The  first  is  near  1.4  -  1.5p  which  could  be  H20,  OH 
first  overtone,  or  CO  second  overtone.  The  second  feature 
appears  near  2.3y  which  could  be  the  CO  first  overtone  observed 
in  Figure  3  for  the  0  ♦  C^  system.  The  third  feature  occurs 
near  2 . 7  y  and  could  be  either  the  band  or  the  OH  fundamental. 

Precise  identification  of  the  emitting  species  from  the  spectra  alone 
is  difficult  because  of  several  reasons.  First  the  system  is  highly  non 
equilibrium  in  nature  with  the  vibrational  temperature  of  the 
emitting  species  far  exceeding  the  rotational  and  kinetic 
temperatures.  Thus  equilibrium  spectra  cannot  be  used  for  easy 
comparison.  Second,  the  overlapping  of  spectra  from  two  or  more 
species  which  is  certainly  expected  in  the  case  of  OH  and  H20 

permit  unique  identification  of  band  shapes.  Third, 
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Figure  3.  The  first  overtone  of  CO  generated  by  reaction  of  0  atoms 
with  in  a  fuel  rich  system  at  1  torr  total  pressure  using  a 

Jarrel  Ash  quarter  meter  monochronator  with  lOOOy  slits. 

1000^  SLITS,  lOO^v  F.S.,  t  =  3,  2/5  RPM  MOTOR,  2  MIN/IN,  71  C2H2, 

16y  0,  808^  He,  200^  02  1st  WINDOW. 
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the  relatively  low  resolution  spectra  do  not  permit  exact 
location  of  individual  line  positions,  therefore,  in  addition 
to  the  spectra  themselves,  we  have  employed  the  kinetic  arguments 
previously  discussed  in  order  to  plan  and  interpret  changes  in 
the  spectra  as  a  function  concentration  nnd  additives. 

First,  we  expect  that  the  spectrum  in  a  large  excess  of  0 
atoms  and  for  small  concentrations  of  HC  (<ly  Malki,  1972) 
should  show  a  predominance  of  OH  emission  over  H20  emission,  at 
least  for  short  contact  times.  In  addition,  CO  emission  should 
also  be  predominant  under  these  conditions  (low  KC,  short  contact 
times)  since  CO  is  a  direct  product  of  the  reaction  fchich  produces  OH 
CHO  +  0  ->•  CO  +  OH. 

It  is  presumably  unlikely  that  CO  is  a  new  bond  in  this  reaction 
and  thus  we  initially  did  not  expect  to  see  the  CO  first  overtone 
from  ethylene  oxygen  atoms  unlike  the  acetylene  plus  oxygen  atom 
case.  Water  is  produced  by  a  subsequent  reaction  of  OH  with 
another  hydrocarbon  fragment.  Hence  H20  emission  should  be 
quadratically  dependent  on  the  hydrocarbon  concentration. 

To  test  these  ideas,  we  examined  the  spectra  from  very  small 
concentrations  (<ly  )  to  nearly  IOQj  of  ethylene.  In  Figure  4, 
we  illustrate  spectra  taken  at  constant  resolution  and  a  constant 
contact  time  of  about  three  milliseconds  over  an  ethylene 
concentration  range  of  2y  to  85y  .  As  illustrated  in  Figure  4, 
the  2 . 7y  peak  predominates  over  the  2.3y  peak  at  low  concentrations 
but  the  reverse  is  true  at  high  concentrations.  In  addition,  there 
are  two  features  at  2.3y  which  also  undergo  an  intensity  reversal. 
Clearly,  a  quadratic  dependence  of  the  intensity  is  not  observed. 


Figure  4.  The  near  infrared  emission  spectrum  generated  by  the 
reaction  of  0  atoms  with  ethylene  for  va’  ious  concentrations  of 
ethylene  at  a  contact  time  of  less  than  2  msec  and  an  oxygen  atom 
concentration  of  approximately  lOy. 

STARTING  FROM  BOTTOM  CURVES: 

2C-H.  (1)  100()  SPLITS,  20,  v  F.S.,  10  F.S.,  x  =  10,  5  min/in, 

m  “  “  y  y  y 

1  RPM,  200y02,  800^He  %  11^0.,  1st  WINDOW 

3  CJ.  (2)  1000,,  SLITS,  20, ,v  F.S.,  10  F.S.,  t  =  10,  5  min/in,  1  RPM 

y  “  y  y  y 

200.,  0?  800,,  He,  Ml  0.,  1st  WINDOW 

y  2  y  y 

10yC2H4  (3)  1000^  SLITS,  50^v  F.S.,  100^  F.S.,  x=10,  5  min/in,  1  RPM 

800^  He  200y02,  11^0,  1st  WINDOW 

45., C9H.  (4)  1000  SLITS,  50  v  P.S.,  100  F.S.,  x=10,  5  min/in,  1  RPM 

y  ^  y  y  y 

800  He  200  0.,  11  0,  1st  WINDOW 

y  y  2  y 

85^C2H4  ,5)  1000y  SLITS,  50^  v  F.S.,  100^  F.S.,  t=10,  5  min/in,  1  RPM 

800  He  200  0, ,  11  0,  1st  WINDOW 

y  y  l  y 
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Figure  4 
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The  feature  at  2.7p  is  very  similar  to  the  observed  OH 
intensity  from  the  reaction 
H  +  03  -  0H++  02 

obtained  by  Murphy  (1971),  illustrated  in  Figure  5.  That 

reaction  is  77  Kcal  exothermic.  The  reaction  of  CHO  +  0  is  even 
more  exothermic  and  hence  many  overlapping  transitions  from  OH 
Av  *  1  are  possible  in  both  cases.  A  cleaner  spectrum  of  OH 
in  which  a  maximum  of  two  vibrational  quanta  can  be  excited 
can  be  obtained  from  the  reaction 
H  +  N02  -  NO  ♦  OH-1* 

At  concentrations  of  N02  sufficiently  low,  the  primary  step  occurs 
in  a  few  milliseconds,  while  the  recombination  of  OH  to  yield 
water  requires  tens  of  milliseconds ,  presumably  via  the  reaction 
0H*+  0H-H20  +  0. 

In  Figures  6  and  7  we  illustrate  the  spectrum  for  0  +  C2H4  and 
the  spectrum  for  H  +  N02  taken  at  (H)  ^  (N02)  ^  4.5.  Two 
overlapping  Av=l  transitions  are  possible  at  2.7u,  but  only 
one  Av»  2  transition,  with  a  resulting  narrower  peak  is  possible  . 
at  1.4  -  l.Sp.  Again  in  this  system,  a  quadratic  dependence  of 
the  intensity  of  H20  emission  on  N02  concentration  would  be 
expected  and  our  kinetic  studies  do  not  show  this.  In  any  event 
no  peak  zt  the  2.3p  region  occurs  in  the  H  ♦  N02  system  and 
therefore  we  conclude  that  the  2.3ppeak  is  no  OH  emission,  but 
rather  CO  first  overtone  emission.  Using  the  NO  4  0  calibration,  we 
find  that  the  ratio  of  the  1 . 5y  intensity  to  the  2.7u  intensity  is 
0.447  compared  to  the  value  of  0.44  *  .03  determined  by  Murphy  (1971) 
for  the  OH  overtone  to  fundamental  intensify  ratio. 
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Figure  5.  The  emission  spectrum  in  the  fun  ^mental  and  first 
overtone  of  OH  due  to  the  reaction  H+Oj  from  Murphy,  (1971). 


Figure  2.  A  comparison  of  the  infrared  emission  spectra  for 
the  reactions  0  ♦  N2I§4 .  0  ♦  C2II4  and  II  ♦  NO,.  Concentrations 
oi  0  and  II  were  approximately  lOu  and  the  concentrations  of 
2**4  *  ^2**4  wcrc  approximately  2u,  lOu,  and  4 .  Su 

rcspcctivrly .  All  spectra  were  taken  at  a  contact  time  of 
less  than  2  msec.  Both  0  and  II  were  diluted  in  a  stream 
consisting  primarily  of  He. 

II  ♦  NO,  -  ;  2000u  slits,  SOuvFS,  lOOuFS,  t  -  10,  1  RPM,  5  Min/i 
SOOullc ,  2  0  0  u  1 1 , ,  SuNO,,  5ul  1st  WINDOW. 

0  ♦  C2II4  -  ;  lOOOu  slits,  SOuvFS,  lOOuFS,  t  -  10,  1  RPM, 

S  Min/ in,  10uC2ll4,  800u!le,  200uO,,  liuO  1st  WINDOW. 

0  ♦  K2II4  -  ;  2000u  slits,  20uvfs,  lOOuFS,  t  •  3,  2  Min/in, 

4  RPM,  2uN2H4  ,  900ullc,  90uO,;  lluO;  1st  WINDOW. 


1.0  2.0  3.0 

WAVELENGTH,/* 

Figure  6. 


Figure  7.  A  comparison  of  the  infrared  emission  spectrum  of 
0  *  *  2^4  ant*  ^  *  NO2  •  Conditions  are  similar  to  those  of 
Figure  6  except  that  pure  oxygen  was  discharged  to  obtain  0 
atoms  and  no  He  was  added.  Hence  a  great  excess  of  0^,  exists 
in  this  system. 

H  ♦  N02  -  :  4.iyK02.  20<M12,  8U0yHc,  4yH2  2000ysl  its ,  100yFS, 
20yvFS,  1  RPM  T  •  10,  1st  K  NDOK. 

0  ♦  C2H4  -  :  2000ysl its ,  630y02,  4.SyC2H4  100yFS, 

100  vFS,  t  •  3,  1  RPM,  1st  WINDOW, 

w  in 
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An  additional  t**st  for  these  tentative  identifications  night 
be  achieved  by  adding  a  reactant  which  preferentially  eliminates 
one  of  the  possible  radiators  by  reaction  or  energy  transfer. 

One  candidate  is  CO,  which  undergoes  a  moderately  fast  reaction 
with  OH  but  is  not  expected  to  alter  H20  emission.  Energy 
transfer  by  excited  CO  to  added  "cool"  CO  a. ght  also  be  expected 
to  alter  the  CO  intensity  distribution.  This  effect  is  readily 
observed  by  addition  of  CO  to  the  oxygen  atom  plus  acetylene 
system  in  Figures,  where  290u  of  CO  is  added  to  the  reaction 
aixturcs.  The  total  intensity  of  the  CO  first  overtone  is 
greatly  reduced  and  the  distribution  is  altered  slightly.  As 
illustrated  in  Figure  9,  an  even  more  pronounced  effect  occurs 
at  the  2.3w  peak  in  the  0  ♦  C2H^  spectrum  at  the  same  CO  concentration. 
At  this  concentration  of  CO,  little  effect  on  the  peaks  at  1.4 
and  2.7u  was  observed.  With  a  rate  constant  of  l.S  X  10  13  c'/scc 
for  the  reaction 

CO  ♦  OH  ♦  CO,  ♦  H 

a 

we  would  expect  an  effective  first  order  rate  constant  of 
l.S  X  10*13  •  l.S  X  103, 

for  the  loss  mechanism  of  OH  via  this  reaction.  This  must 
be  compared  to  the  oth<  r  loss  mechanism  for  OH+  , 

OHf  ♦  0-02  *  H, 

which  has  a  rate  constant  of  approximately  4  X  10  11  cc/scc.  Hence, 
with  an  oxygen  atom  cor  cent  rat  ion  of  approximately  3  X  10l4/cc 
the  dominant  OH  loss  mechanism  remains  0  atom  attack,  with  an 

4 

effective  first  order  rate  constant  of  1.2  X  10 

In  the  ease  of  the  H  ♦  N02  system,  where  OHf  loss  is 


figure  8.  The  effect  of  an  addition  of  150p  of  CO  to  the 
0  ♦  C2H,  system. 


Conditions  for  Both;  1000yslits,  £  RPM,  -5- -t^-n ,  T  *  10, 
SOyVFS  S9yC2H2  ,  208y02,  812  He,  ^8..0  270..  Ar. 


3 


Figure  9.  The  effect  of  additives  on  the  0  +  C2H4  emission 
spectrum.  240^C09  and  236^CO  were  added. 

Conditions  for  all:  1000  slits,  1  RPM,  5>Min,  t  =  10, 

y  ’  in  ’  ’ 

50  vFS,  100  FS,  50  C-H .  ,  208  0,  812, He,  8,0,  1st  (KB  )  WINDOW. 
M  y  p  “  ^  y  l  y  y  r 


1.0 


2.0 

WAVELENGTH,  A 


3.0 


l 

I 


Figure  9. 


presumably  via 


0H+  +  OH  -*•  H20  +  0, 

a  much  slower  reaction,  it  should  be  possible  to  attenuate  the 
OH  signal  via  CO  scavenging.  As  shown  in  Figure  10  this  indeed 
does  happen.  In  addition,  other  species  which  react  with  OH 
have  been  added  to  this  system.  In  general,  we  see  from  Figure  10 
that  the  more  rapid  the  reaction  with  OH,  the  smaller  the  signal 
in  both  the  OH  fundamental  and  the  OH  overtone  .  The  overtone 
decreases  more  rapidly  than  the  fundamental,  which  is  con¬ 
sistent  with  the  narrowing  of  the  fundamental  band  also  observed. 

Thus,  we  conclude  on  the  basis  of  kinetic  arguments  as 
well  as  spectral  arguments  that  the  observed  emission  at 
2.7y  is  due  to  the  OH  fundamental  and  the  emission  at  1.5y 
is  the  OH  overtone. 

As  shown  in  Figure  11,  we  also  observe  these  features  in 
the  H  +  HO 2  system;  however,  the  conditions  necessary  to  achieve 
this  signal  level  in  Figure  11  make  it  impossible  to  permit 
H  +  H02  to  be  the  cause  of  the  OH  emission  in  the  0  +  C2H4 
system. 

In  summary,  we  present  in  Figure  12,  a  compilation  of 
infrared  signatures  at  a  contact  time  of  2  msec  due  to  0  atom 
attack  on  various  hydrocarbons. 


V 


Figure  10.  The  effect  of  additives  on  the  H+N02  infrared  emission 
signal . 

Conditions  Common  to  All:  2000^slits,  ,  1  RPM,  774^ile,  222^ 

H_ ,  5  NO?  1st  WINDOW,  t=10,  100  vF.S. 
z  y  z  m 

Starting  from  Top: 

N02+H  -*•  (Ar):  1  Torr  Ar ,  1  Torr  F.S. 

NO  +H  -*•  (CO):  780  CO,  1  Torr  F.S. 

z  y 

N02+H  ■>  (C2H2):  310^C2H2,  1  Torr  F.S. 

N02+H  -  (C2H6):  0.5  Torr  C2H6,  1  Torr  F.S. 

N02+H  -  (n-C4H1Q):  370yn-C4H1Q,  1  Torr  F.S. 

N02+H  -  (C3H6):  345yC3H6,  1  Torr  F.S. 
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Figure  11.  A  comparison  between  the  H+N02  and  H+HO,  infrared 
emission  signals.  Exactly  320Gy  02  was  added  to  the  in  He 
discharge  in  order  to  observe  the  H+H02  signal  at  a  contact  time 
of  2  msec. 


Conditions 
F.S.,  IRPM, 


for  Both: 

5  min 
in  ’ 


2000  Slits,  T-10,  1st  WINDOW  100  F.S.,  100  v 

M  p  p 

222  H,,  774  He.,  5  NO,. 

V  ‘  P  m2 
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i-igurc  12.  A  Suwary  of  infrared  signatures  initiated  by  0  atom 
attack  on  selected  hydrocarbons. 


Conditions  for  All:  1000uSlits,  IRPM  t-10  SO 

10°UF.S.,  1st  KINDOK,  774uHe  195^02,  250yAr. 

Starting  fro»  Top: 


Cjll  *0 

4  • 

• 

c2i.4.° 

-*  : 

S5PC2H4 

Cj'V° 

-*  : 

70uC3H6 

n-Cn"lO 

♦0  •* 

:  ‘VC4H 

c2h6.o 

4 

• 

36UCA 
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THEORETICAL  DETERMINATION  OF  VIBRATIONAL  EXCITATION  OF  MOLECULES 

The  theoretical  program  at  RIES  has  been  primarily  concerned 
with  the  calculation  of  cross  sections  for  the  vibrational  ex¬ 
citation  of  molecules  by  heavy  body  collisions.  This  is  carried 
out  using  a  numerical  close  coupling  solution  of  the  wave 
equation  for  the  colliding  system  based  on  the  breathing  sphere 

i 

model  of  Schwartz,  Slawsky  and  Herzfeld  (1954)  to  represent  the 
target  molecular  vibrator,  and  an  empirical  spherically  symmetric 
scattering  potential. 

M?ny  reactions  of  interest  in  the  upper  atmosphere  involve 
collisions  between  molecules  and  atomic  oxygen.  However,  except 
for  a  semi-empirical  estimate  of  the  CO-O^D)  interaction 
(Marriott,  1970)  the  central  potentials  necessary  for  the 
calculation  of  these  excitation  cross  sections  have  not  been 
available  until  recently. 

As  detailed  in  Part  III  of  this  report,  potentials  of  the 

form 

V(r)  *  Aexp(-Br) 

have  now  been  determined  experimentally  for  a  number  of  col¬ 
lision  systems  of  interest.  The  parametric  values  required  by 
the  collision  code  are  shown  in  Table  1. 

Up  to  the  present  time,  using  this  data  excitation  cross 
sections  have  been  calculated  for  the  C0-0(3P)  and  N2-0(3P) 
systems.  Preliminary  results  based  on  the  coupling  of  three 
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molecular  vibrational  states  are  shown  in  Table  2,  together  with 
the  CO-CO  and  CO-O^D)  values  for  comparison. 

Although  no  direct  cross  section  measurements  are  available 
for  checking  these  theoretical  results,  experimental  determi¬ 
nations  of  certain  rate  coefficients  are  available  (center- 
1973)  and  these  can  in  turn  be  calculated  using  the  cross 
sections.  Thus  an  indirect  comparison  with  experiment  can  be 
obtained  and  this  calculation  is  presently  underway. 

In  addition  to  the  above,  the  previous  report  outlined 
a  recent  extension  of  the  original  vibration- translation  col¬ 
lision  code  to  include  vibration-vibration  exchange  processes. 
The  need  for  further  extension  to  include  virtual  states  was 
established  if  the  code  was  to  continue  to  be  able  to  calculate 
cross  sections  at  collision  energies  of  interest  (i.e..  below 
1-2  eV) . 

The  code  has  now  been  reconstructed  to  allow  for  the 
occurrence  of  these  virtual  states.  At  this  time  the  new 
subroutines  are  being  debugged. 

Our  present  plans  for  future  work  include  the  following: 

(a)  Complete  the  debugging  of  the  new  v-v  code 
including  the  effects  of  virtual  states. 

(b)  Complete  the  vibrational  excitation  cross 
section  calculations  for  a?  the  collision 
systems  for  which  0-molecule  scattering 
potentials  are  now  available  (Table  1). 


(c)  Apply  the  full  code  initially  to  a  detailed  study 
of  vibrationa.  excitation  rpocesses  in  CO,  in¬ 
cluding  the  effects  of  vibration-vibration  exchange 
and  virtual  states.  The  results  will  be  analysed  in 
comparison  with  previous  theory  and  experiment  to 
estimate  the  "physical"  reliability  of  the  overall 
calculation . 


Table  1.  Experimental  central  potential  parameters  for  the 
interaction  V(r)  =  Aexp(-Br)  for  various 
0(JP) -molecule  collisions. 


System 

A  (atomic  units^ 

B  (ao1} 

Expt.  range 

O  -  CO 

30.77 

1.7447 

13.6  -  0.4 

0  -  N2 

20.92 

1.6108 

14.8  -  0.4 

0  -  NO 

38.13 

1.7717 

13.4  -  0.4 

0  -  H2 

18.38 

2.1474 

14.9  -  0.4 

0  -  02 

62.43 

1.9151 

12.4  -  0.4 

0  -  C02 

75.55 

1.7214 

12.4  -  0.4 

0  -  h20 

52.02 

1.9579 

12.9  -  0.4 

Table  2.  Calculated  cross  sections  for  the  vibrational 

excitation  of  N„  and  CO  assuming  coupling  between 
the  first  three  vibrational  levels. 


Collision  system: 

Nn  -  0(3P) 

4  2 

CO  -  0(3P) 

CO  -  co+ 

CO  -  0(1D) 

Coll,  energy  (eV) 

Q0i 

Q01(lTa0) 

Q01(ira0} 

Q01("a0> 

0.266 

0.292 

0 .000 

0.000 

0.000 

0.000 

0.6 

0.117 (+2) 

0.8 

0 . 578  (-8) 

0.410 (-6) 

0 . 680 (-5) 
0.339 (-3) 
0.265 (-2) 

0.85 

1.0 

0. 178  (-5) 
0. 214  (-4) 
0. 277  (-3) 

0 . 56 (-4) 

0 . 110 (+2) 
0 . 101 (+2) 
0.85  (+1) 
0.70  (+1) 

1.5 

2.0 

0 . 12  (-1) 

2.5 

0.22 

0.60  (+1) 

**mistrbeiLS?pned.Sh°WS  P°“er  °f  ten  by  Which  number 


tMarriott  §  Gianturco  (1970) 
ttMarriott  (1970) 
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PART  III:  Determination  of  the  Intermolecular  Potentials 
between  Oxygen  Atoms  and  Plume  Species. 


P.K.  Rol 


A  molecular  beam  apparatus  donated  by  NASA  has  been  used 
to  measure  intermolecular  potentials  required  as  input  for 
close  coupled  calculations  of  excitation  cross  sections. 

Apparent  cross  sections  for  different  energies  are 
determined  by  measuring  the  attenuations  of  a  fast  (50  -  2500  eV) 
atom  beam  in  passing  through  a  scattering  gas.  The  interaction 
potential  is  calculated  from  these  cross  sections  and  the 
scattering  geometry. 

The  molecular  beam  apparatus  used  in  this  experiment  is 
shown  in  Fig.  iva.  Atomic  ions  are  produced  in  the  ion 

source.  The  ion  beam  is  analyzed  and  the  selected  ions  are 
forcused,  collimated  and  decellerated  to  the  required  energy 
prior  to  neutralization  in  a  charge  exchange  chamber.  Just 
before  entering  the  charge  exchange  chamber  the  ion  beam  is 
deflected  to  prevent  high  energy  neutrals,  formed  by  neutra¬ 
lization  with  the  background  gas,  to  be  mixed  with  lower  energetic 
atoms.  The  neutralized  beam  is  finally  passed  through  the 
scattering  chamber  and  detected  with  a  channel  tron  multiplier. 

The  beam  intensity  is  measured  as  a  function  of  the  scattering 

gas  pressure,  and  the  cross  sections  S  can  be  calculated  ac¬ 
cording  to: 


Where  I  =  attenuated  beam  flux. 

I0  =  incident  beam  flux, 
n  =  scattering  gas  density, 
n  l  =  path  length  in  scattering  gas. 
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p  — 


The  equipment  was  tested  by  measuring  the  He -He  potential  for 
which  accurate  theoretical  calculations  are  available.  Ex¬ 
cellent  agreement  was  obtained  with  these  theoretical  results. 

The  repulsive  potentials  for  oxygen  atoms  and  several  atoms 
and  molecules  have  been  calculated  from  measured  incomplete  scat¬ 
tering  cross  sections  of  fast-atom  beams.  These  fast  beams  were 
produced  by  neutralizing  either  positive  or  negative  0-ion  beans 

m  a  charge  exchange  gas.  H20,  Xe,  02  and  Kr  were  used  as  charge 
exchange  gases. 

The  results  are  given  in  the  lorm  of  Born-Mayer  potentials 
in  the  following  table. 

Oxygen-atom  Potentials 

-br  I  A  in  eV 

Ae  r,<r<r  .  o.i 

1  u  jb  in  X  1 

(  rl -  ru  in  X 


Species 

A 

b 

rl 

ru 

Exchange 

Reactions 

N2 

569 

r 

3.044 

1.20 

2.37 

o\  h7o 

CO 

837 

3.297 

1.25 

2.32 

o  ,  h7o 

C°2 

2055 

3.253 

1 . 57 

2.64 

o  ,  h7o 

H2 

500 

4.058 

0.865 

1.75 

0  ,  Xe 

NO 

1037 

3.348 

1.30 

2.35 

0+,  Xe 

h2o 

1415 

3.700 

1.27 

2.22 

0+,  Xe 

°2 

1698 

1  3.619 

l .  36 

2.32 

0+,  Xe 

He 

355 

3.728 

0.833 

1.80 

o',  o? 

Ne 

1415 

3.984 

1.18 

2.06 

0- ,  0 

At 

745 

3.184 

1.25 

2.36 

2 

o',  o7 

Kr 

1040 

3.239 

1.34 

2.43 

O',  0 

Xe 

1470 

3.303 

1.44 

|  2.50 

o',  02 

He 

360 

3.702 

0.844 

1.82 

0+,  Xe 

Ne 

1312 

3.988 

l.lu 

2.04 

0+,  Xe 

Ne 

1125 

3.864 

1.15 

2.06 

0%  Kr 
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m  addition  the  incomplete  cross  sections  have 

the  following  combinations  and  potentials 
eulated: 


been  measured 
ore  being  cal- 


Hc-Hc 

He-Ne 

He-Ar 

He-Kr 

Hc-Xe 

and  H-Hc,  D-Hc. 


Ne-Ne  Ar-Ar 

N’c-Ar  Ar-Kr 

Nc‘Kr  Ar-Xe 

N'e-Xc 


Xc-Xc 

Kr-Xc 


52 


Figure  IVa.  Amdur-type  molecular  beam  apparatus  at  Wayne  State. 


